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^ ! Abstract 
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f^ . A renormalization scheme for the nucleon-nucleon (NN) interaction based on a sub- 

CN I tracted T-matrix equation is proposed and apphed to the one-pion-exchange poten- 

tial supplemented by contact interactions. The singlet and triplet scattering lengths 
Q^ ■ are given to fix the renormalized strengths of the contact interactions. With only 

one scaling parameter (/i), the results show an overall very good agreement with 
neutron-proton data, particularly for the observables related to the triplet channel. 
The agreement is qualitative in the ^Sq channel. Between the low-energy NN ob- 



^ , servables we have examined, the mixing parameter of the "^Si — "^ Di states is the 
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most sensible to the scale. The scheme is renormalization group invariant for /i -^ cxo. 
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1 Introduction 



The effective field theory (EFT) of nuclear forces based on a chiral expansion of an effective 
Lagrangian, as proposed by Weinberg [1], gives a nucleon-nucleon (NN) interaction in 
leading order which is the one-pion-exchange potential (OPEP) plus a Dirac-delta. He 
suggested to infer the values of the strength of the delta force in the ^5*0 and ^5*1 channels 
from the singlet and triplet scattering lengths respectively. Also such potential should be 
valid for momenta well below a typical momentum scale of quantum chromo dynamics such 
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as the rho meson mass {rrip ~ 4 f m ^) [1]. Therefore, the intermediate virtual propagation 
of the NN system should be cut at the momentum scale of this order. 

In spite of the significative results obtained by several authors, when implementing the 
EFT program for the two nucleon system [2-8], and also by a recent related approach [9], 
as far as we know the predictibility power of the leading order term (OPEP plus delta) 
with a single renormalization momentum scale was not fully explored in the literature. 
The result of such calculation would give a sound basis for the renormalization program 
of EFT in NN system. The leading order term was already considered in NN calculations, 
when comparing with more realistic calculations [10], or as part of interactions used in 
the EFT program [3,4,11,12]. However, we argue that, without a systematic analysis of 
the physical contribution coming from each order term in the renormalization procedure, 
it is easy to miss their real significance when parametrizing the full interaction to obtain 
the desired observables. To start, it is relevant to know how much physical information 
can one extract from the leading order term, using a renormalization procedure. What are 
the observables one can better describe with just the leading order term (the higher order 
terms should not affect too much these observables) and what are the observables for which 
higher order terms are essential? Our aim in this paper is to answer such questions. We 
hope the following analysis can be useful in the renormalization program of the effective 
field theory applied to the nucleon-nucleon system. 

In this work we have verified strictly the validity of Weinberg suggestion, to first consider 
(in leading order) an effective potential which "consists of just a conventional one-pion 
exchange term, plus a direct neutron-proton (n — p) interaction produced by the four- 
fermion terms" [1]. Such "four- fermion terms" consist of two delta- function interactions: 
one in the ^5*0 channel and the other in the ^5*1 channel. To be consistent with the 
original suggestion, to cut the virtual propagation of the NN system at the scale rup, 
we first use a renormalization scale of such order. However, we understand this as a 
qualitative suggestion, such that we explore its limitation in the calculation of several 
low-energy NN observables. So, we finally consider a free renormalization scale in our 
numerical calculation, and try to obtain a consistent description, in order that such scale 
vanishes as a physical parameter as required by a renormalization group invariant theory. 

In our approach, we use a renormalization scheme for the T-matrix, which is appropriate 
when considering the OPEP plus a delta potential. It incorporates the momentum scale 
(/i) where the two-nucleon propagation is removed and the scattering lenghts fix the 
renormalized strengths. A subtracted T-matrix integral equation is derived, where the 
subtraction is done at energy of — /i^, in the center of mass system (in units oi h = 1 
and nucleon mass m = 1). The input is the T-matrix at energy — /i^, which is chosen 
as the OPEP supplemented by delta interactions. The scale fi is of the order of rrip, and 
we will show that the results for the renormalized T-matrix, in the limit of /i — » oo, once 
the scattering lengths are kept fixed, are quite consistent with /i ~ nip, at least for the 
diagonal components of the matrix (the extension of this statement will be clarified). As 
we verified, the numerical results in the ^5*1 — ^ Di channels, obtained by following this 
procedure are in good qualitative agreement with the experimental data for the n — p 



phase-shifts, mixing parameter, the deuteron binding energy and D/S ratio, at laboratory 
energies below ~ 50 MeV. 

The higher order terms in the EFT expansion of the effective interaction will become more 
important for laboratory energies above 50 MeV in the ^5*1 —^ Di channels, where the 
leading order results clearly deviate from the data. However, in the ^5*0 state the impor- 
tance of the higher order terms appears already in the value of the effective range, which 
is underestimated in lowest order. Such results confirm that the effective interaction in 
the singlet state needs higher order correction (as already expected), but such corrections 
should not affect in an essential way the triplet states (particularly at low energies), which 
are already in a fair agreement with data. 

In our scheme, it is considered an unregulated effective power expansion of the NN poten- 
tial only in the mid- and short- range parts of the interaction, such that the well stablished 
long-range part (OPEP) is kept intact Q. So, the matrix elements of the full effective NN 
interaction, in the relative momentum space, will be given by 



{p'\Veft\p) = {p'\V^\^ + 
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where the unregulated strengths of the contact interaction, and the higher order deriva- 
tives, have subindices t and s, for the isospin singlet (spin triplet) and isospin triplet (spin 
singlet) channels respectively. The matrix element of the one-pion-exchange potential is 
given by: 



{p'\vM 



9l 



4(27r)3/; 



:Ti ■T2- 



<?1 ■ JP' -P) ^2- JP' - P) 



(2) 



where ai and Tj are the usual spin and isosping Pauli matrices for nucleon z; Qa (= 1.25) 
is the axial coupling constant, /,r (= 93 MeV) is the pion weak-decay constant, and 
ttItt (= 138 MeV) is the pion mass, ff — p is the momentum transfer. The leading order 
term of the expansion, which we are going to consider in the next sections, imply to take 
A^?A = for all n > 1 in eq. (1). 

In section II we present the formalism for the subtracted T-matrix, which is applied to 
the leading order term of the NN interaction. In section III we present our main results 
and conclusions. 



^ The OPEP is common to all realistic interactions, and an effective power expansion which 
does not consider it will have to face difficulties with the parametrization, as the constants 
should represent simultaneously long and short range physics. 



2 Subtracted T-matrix and the leading order NN interaction 



In the next we introduce our formalism which is based on a subtracted equation for the 
T-matrix. We begin with the formal T-matrix equation (for a regular potential), which is 
written as: 

T{E) =V + VG\^\e)T{E) = [1 - VG\^\e)]-W (3) 

= V + T{E)G^^\e)V=V\1-G^^\e)V]-^ , (4) 

where the free Green's function for the two-body system, with the appropriate boundary 
condition, in terms of the free Hamiltonian Hq, is 

Gi;\E) = {E + te-Hor\ (5) 

Considering that our input is given by the T-matrix at a given scale, the potential will 
formally be replaced by the T-matrix at scale n through Eq.(3) or Eq.(4): 

V = T{-f,') [l + Go( V)n V)] ~' = [l + T{-f,')Goi-fi')] ' n-f,') . (6) 

Substituting the potential given by Eq.(6) in Eqs. (3) and (4), the corresponding sub- 
tracted equation for the T-matrix is found: 

TiE) = T( V) + T( V) {Gi''\E) - Goi-f^')) TiE) (7) 

One should note that we have chosen a real value for fi (negative energy), but this is 
not a restriction to the generality of the subtracted T-matrix equation (7). However, such 
choice has the advantage that the kernel of the integral equation has just one fixed point 
singularity for positive energies E, and the same formal structure as the original T-matrix 
equation. As one can see, from the original T-matrix, one simply needs to replace V by 
T{—fi^) and multiply the free propagator by a function which depends on the scale n, 
such that 

T{E)=Ti-^^')+Ti-^^')G^^\E■,-^^')TiE) (8) 

where 

gW(E; -^^) = Gi^\E) - G,{-^?) = l^^fi^\E). (9) 

The renormalized T-matrix equation, given above, has a kernel which is subtracted at 
the energy point — /i^. Although the interaction V, as given by Eq. (6), is not defined 



for a singular interaction, the subtracted T-matrix equation is well defined for a regular 
interaction plus a delta potential. This justifies a posteriori the formal manipulations in 
obtaing Eq.(7). One could instead, work with a finite cut-off in the momentum integrations 
and derive the subtracted equation for the T-matrix; then take the infinite limit of the 
cut-off. Of course the subtracted equation for the T-matrix is insensible to the cut-off as 
long the potential is singular like a delta interaction. 

In the second term of the right-hand side of Eq.(7), the propagation over the intermediate 
states is subtracted at an energy scale of — /i^. The physical information apparently lost 
in this subtraction is contained in T(— yU^), which in the case of pure contact interaction 
is the renormalized coupling constant [13]. 

If one sets T{—fi^) equal to a given operator (/i independent), the scheme becomes depen- 
dent on the subtraction point, which only has /i vanishing as a physical parameter (i.e., the 
theory is renormalization group invariant) in the limit fi -^ oo. This is a trivial observa- 
tion when such operator is a regular potential. Thus, setting the potential as equal to the 
T-matrix at ii^ = — /x^, the present scheme implies, that the physical T-matrix is found in 
the limit ;U ^ oo. We will show that the procedure can be applied for a singular potential 
like OPEP plus contact term, and the limit /x ^ oo gives finite physical observables with 
vanishing dependence on /x. The numerical results obtained by constraining T(— yU^) to be 
equal to OPEP plus contact term for several values of n show quantitatively the limiting 
process, while the scattering lengths are kept fixed. In this way, one can also verify how 
good is a scale as low as the p-meson mass, in order to approach the limiting values for 
the observables. Anticipating our numerical findings, the low-energy phase-shifts for the 
^5*0, ^5*1 and ^Di channels give a reasonable support to such low value of fi. However, 
this is not true for the ei mixing-parameter, even for energies as low as 20 MeV. 

Let us apply the subtracted T-matrix equation for the pure contact interaction, and re- 
obtain some well-known results just for the sake of completeness. In this case, the matrix 
elements of the T-matrix are independent of the momentum and, at E = — /i^, it is given 
by the renormalized coupling constant: 



(p'|T(V)b-) = ^A^(/^). (10) 

The renormalized T-matrix in the subtraction point is given by X-jilfi), which is the 
physical information supplied by the two-body system in order to obtain other observables. 
For example, A7^(0) is identified with the scattering length (a), A7^(0) = a. Solving Eq.(7), 
in the usual three-dimensional space with condition (10), 



{p'\Tmp) 



The on-shell-momentum is k 
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{p'\T{Em 
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A7^(/i) 



+ (/i + ik) 



:i2) 



In general the matrix elements of a two-body interaction will depend on the momentum, 
spin and isospin, such that it will be usefull to define our notation for the partial wave 
decomposition. The decomposition of the plane-wave is given by 



V TT , . r 



(13) 



where the quantum numbers j,mj,l,s,I refer, respectively, to the two-body operators: 
total angular momentum J, its z-component Jz, orbital angular momentum L, total spin 
S and total isospin I. yl^^ip) are the usual orthonormalized functions for the coupled 
angular (L) and spin (S) momentum to the total momentum (J = L + S). p defines the 
two-dimensional angular variables related to the vector p. 

The partial wave decomposition of OPEP (2) and, for simplicity, dropping the unnecessary 
indices for the total angular momentum and spin in VJr, is given by 



V^'I^ij^.p) 



,^[2/.(/. + l)-3]^^.. .,/^, xt^ 



(Ji ■ {p' - p)a2 ■ {p' -p) 



y;i'{p'] 



{p' — p)"^ + ml 
where for the singlet state a = s {Ig = 1) and for triplet state a = t {It = 0). 



(14) 



So, in the angular projected equation for the T-matrix, the singlet ^5*0 state matrix element 
for the OPEP, corresponding to eq. (2), is given as 



vZ'Kp'^p) 



al 



m: 



2 — dx 

327r/^ ^ J p'^ + p'"^ — 2pp'x + ml 



(15) 



In the limit of p or p' going to infinity the matrix element goes to a constant: 



Jim V;(f(p',p) = -f^ 



(16) 



The corresponding angular momentum projected matrix elements, for triplet ^5*1 — ^ Di 

(I'D 

coupled channels, (V^ j , with /', / = 0, 2) are: 



VT iP', P) = W^A^- f dx . ,, 7 ^ ^ , , (17) 

,.(20). , s glV^ }, p'^-2pp'x + p\lx^-\) 

K,* '{.P,P) = T^-F^ / dx ^—^2 o '\ 2 ' 

\/(f (p' p) = ^^l d /PP'^-(p'+p"n^'--2) ng) 



3271 f!^ J p^ + p' — 2pp'x + m^ 

and<f(p',p) = \/if(p,p'). 

In the limit of p or p' going to infinity, the matrix elements of OPEP in the ^^i — ^ Di 
channels are zero or constant: 

lim Vlf (p', p) = -f- , lim K? ) (p', p) = 0, 

^lim V^f iP',P) = §^, Jim <f ) (p', p) = ; (20) 

consequently, just one subtraction is enough to obtain a finite T-matrix, once T(— /i^) is 
given. 

The trivial choice for T(— /x^) is the OPEP operator added to a delta interaction. For the 
singlet renormalized T-matrix, with the renormalized strength of the contact interaction 
{\n,s{fJ')) fixed by the singlet scattering length a^ = —23.7 fm, we have 

T^''\p',p; -^') = Vj'^\p',p) + A^,,(/i) , (21) 

where we have dropped the upper index in A for simplicity. As for the coupled ^Si — '^ Di 
channels, the renormalized T-matrix for each channel (/, /' = 0, 2) is given by 

Tp\p',p; V) = VJi;'\p',p) + Xn,tWi'oSio ■ (22) 

In this case, the renormalized strength of the contact interaction {Xn,t{fJ')) is fixed by the 
triplet scattering length, at = 5.4 fm. 

Next, the formalism in momentum space, for the renormalized partial- wave T-matrix, 
using one subtraction in the propagator at an energy — /i^ is concluded by writting the 
integral equations corresponding to Eq. (8) where one should replace the matrix elements 
for T(— /i^) by the corresponding above expressions [for the singlet state, eq. (21); for the 
triplet state, eq. (22)]. For the singlet state, we have 



Tr\p',P;k') = Tr\p',P;-^') + 



and for the triplet state, with /j = 0, 2 (? = 1,2,3), the corresponding coupled equation 
is given by: 



t('^'^Hp',P;P)=T/'^'^)(p',P;V) + 



2y-7^ 2(f^' + k'\T^'''''\p',q;-n')^i,,,) 



-J2 dqq 

71 , J 



I v/^^ + Q'^/ k?' — q^ + ie 



Tr'>{p',p-e) (24) 



These eqs. (23) and (24) for a regular potential clearly converge to the exact original 
T- matrix in the limit /i^ — > oo, with the corresponding Born terms given by the potential. 



3 Numerical results and conclusions 



In the following, we present the numerical results for the observables, in the ^5*0 and 
^5*1 —"^Di channels, calculated from the leading order NN interaction, with the subtracted 
eq. (7) and the T-matrix elements at scale fi defined by eqs. (21) and (22). Distinct 
numerical procedures can be used to obtain the final results, from the equations (23) and 
(24). In our numerical calculation, we have followed the one described in Ref. ([14]), where 
the T-matrix is related to an auxiliary real matrix, which is solved by the inverse matrix 
approach. Considered that in the present work we are not interested in high numerical 
accuracy, any of the usual numerical procedures should give about the same results we 
obtain, once a large enough number of points is taken in the numerical discretization of 
the momentum space. In our calculation, when discretizing the momentum space (p) from 
to cxD, we have used a Gaussian mesh with up to 100 points, mapping the interval [—1,1] 
to [0, oo] through the expression p = po{l + x)/{l — x). po is adjusted to optimize the 
convergence of the results, and in practice should be comparable with /i, when a minimum 
number of mesh points is required. 

In Table I, for several values of fi, we show the corresponding results for the singlet and 
triplet effective ranges, and also for deuteron observables (binding energies and ratio r^^). 
The ratio rjr) is given by 

TjD = lim —TTT^. , (25) 

^-**=sTf°^(A;,A;;P) 

where the deuteron binding energy is E = —k\. The effective ranges, for the singlet state 
(ro,s) and for the triplet state (ro,t) and the corresponding scattering lengths (a^ and at) 
are related to the s— wave phase shifts, according to the effective range expansion: 



kcot (6i;:^{k)) = -- + Ko^^e + 0{k') . (26) 



Table I Low energy n-p and deuteron observables compared with data. Singlet (rg^s) and 
triplet (ro^t) effective ranges, deuteron binding energy {Bd) and ratio r//) are given for several 
values of the single parameter ^ which is used in our renormalization scheme. The \tz^s{ij) and 
\'R.,t{p) are the strenghts of the 5— interactions which were added to the OPEP and adjusted to 
the corresponding scattering lengths, a^ =-23.7 fm and at = 5.4 fm. The results were compared 
with data obtained from Ref. [15]. 



/i (fm-i) A7^,,(fm) ro,,(fm) A^,i(fm) ro,t(fm) Bd (MeV) t^d 



4 


-0.8806 


1.332 


-0.2281 


1.364 


1.977 


0.02808 


10 


-0.7570 


1.345 


21.741 


1.536 


2.084 


0.02904 


30 


-0.6977 


1.347 


-0.3776 


1.582 


2.114 


0.02933 


Ref. [15 


- 


2.73 


- 


1.75 


2.2246 


0.0256 



The results in Table I are compared with experimental results obtained from Ref. [15], 
and show a good qualitative agreement, considering that we are only using the leading 
order term of the effective NN interaction. As it is shown, the agreement improves as 
we increase the value of /x. The results in the ^5*0 channel could be improved by adding 
higher-order terms in the effective interaction, as already one can deduce from previous 
approaches used by other authors in EFT [4,7,8]. It is also important to observe that for 
/i > 4fm~^, tqs is pratically stable, such that it is insensible to nucleon propagation for 
higher momenta once the scattering length is fixed. The insensibility to /i > 4fm~^ is also 
found in the singlet phase-shifts up to laboratory energies of about 100 MeV. 

In the triplet channel, we observe a large variation of \n,t{lJ') as the scale is moved to higher 
values. The oscillating behaviour is a consequence of the zeroes of X-jz^tif^) at those scales 
where only OPEP is enough to reproduce the triplet scattering length. However, in spite 
of such variation, the triplet low-energy and deuteron observables are quite stable, and 
converging in the limit of /i going to infinity. From the observables we have examined, 
as we pointed out in our discussion, the slowest rate of convergence is found for the 
mixing parameter ei, which is defined through the on-shell S— and T— matrices in the 
coupled channel. With an implicit energy dependence, such relation can be written as a 
two-dimensional matrix and parametrized according to Ref. [16]: 



S = l- 2ikTt =1 ^ ' ^ ' I . (27) 

i sin(2ei)e*(''"'*+''2,t) cos(2ei)e2*'52.t 



In Figs. 1, 2 and 3 we show the results for the neutron-proton phase-shifts ^5*0, ^^i and 
^Di, respectively. The results shown in Fig. 1 are pratically insensible to nucleon-nucleon 
propagation for momenta higher than 4 fm~^, in agreement with results shown for the 
singlet effective range (Table I), such that no improvement can be expected for the singlet 
observables with just one subtraction/parameter. In Fig. 2, the good qualitative fit to 
data obtained with the scaling parameter /x = 4 fm ~^ is improved when taking larger 
value for /i. The subtraction point above /x > 4fm~ has some effect, as also anticipated 
from the values of scattering length and effective range (see Table I); but the effect is 
not strong above 30 fm~^. In Fig. 3, we observe that the insensibility to nucleon-nucleon 
propagation at higher momenta is also found in the D— wave. In this state the phase-shift 
is clearly dominated by the OPEP, as also shown in ref. [10], in a first order perturbative 
calculation. 

In Fig. 4 we present our results for the mixing parameter, ei, which is obtained using 
eq. (27). The results are insentive to /^ > 4fm^ only for Eiab < 10 MeV, and show a fair 
reproduction of the Nijmegen values. This follows due to the dominance of OPEP in the 
slope of ei at zero energy [17]. Above 10 MeV, the mixing parameter shows a dependence 
on the values of /i, which is reduced by increasing /i above 30fm~^, while the agreement 
with data is improved as, ^ ^ oo. The mixing parameter is the most sensitive parameter 
in respect to changes in the subtraction point for Eiab > 10 MeV and for fi between 4 
and 30 fm~^. Our results in the limit of /i ^ oo are consistent with the ones obtained by 
Ballot and Robillota [10], where they have used a cut-off parameter. 

The results presented in Table I and Figs. 1-4, when compared with fitted data [15] 
and with other calculations with more fitting parameters [3,4,11], are a clear evidence 
of the dominance of the leading order term (OPEP) in low-energy NN observables. This 
implies that one should first renormalize the T-matrix for the leading order term and 
then add higher order terms to the NN interaction in the chiral expansion, or in the 
EFT expansion. Our results for the phase-shifts and mixing parameter show that the 
mixing parameter is the most sensitive observable to the subtraction point used in the 
renormalization procedure. The phase-shifts for ^5*0 state and for the ^Di are pratically 
insensitive to /i larger than 4 fm~^, with an evidence that ^Di is dominated by the OPEP, 
in agreement with previous calculations [10]. The present results confirm that the addition 
of higher order singular terms in the potential are needed in order to fit the singlet effective 
range and corresponding phase-shifts. This observation was already verified in several EFT 
approaches [3-7,11,12], without however addressing to the full predict ibility of the leading 
order term. This lead some authors to affect negatively the nice fitting obtained for the 
triplet channel when trying to improve the observables for the singlet channel [4] . In order 
to avoid this problem, it is important to parametrize the higher order corrections, such 
that it affects mostly the singlet channel (where the leading order fails to give a good 
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description also to the effective range, independently of the renormalization point). In the 
^5"! —^ Di channels, the higher order terms will become more important for laboratory 
energies above 50 MeV, where the leading order results deviate from the data (as shown 
in Figs. 2 and 4). 

Considering the results obtained and shown in Figs. 1-4, we would like to make some 
comments about the renormalization scale /i, which appears in our approach as the point 
where we introduce a subtraction in the nucleon-nucleon propagator. (This parameter is 
directly related to the "cut-ofF' parameter of other approaches.) Weinberg [1] has origi- 
nally suggested to cut the virtual propagation of the NN system at the scale nip, which 
implies /i ~ 4 fm~^, and this suggestion is consistent with our results, as it is almost 
exactly verified for the ^5*0 and ^Di phase shifts (see Figs. 1 and 3), and also approxi- 
mately for the ^5*1 phase shifts (Fig. 2). We further observe that our results for fi —>■ oo 
(numerically, we have used /x < 60 fm~^) are about the same we got for the ^5*0 and ^Di, 
and improves the results we obtain for the ^5*1 case, for the mixing parameter, and also 
related observables, as shown in Table I. Such results, particularly the results for the mix- 
ing parameter, lead us to conclude that Weinberg's suggestion can only be understood as 
a qualitative statement, and as such, it is verified. However, we observe that the mixing 
parameter, which is very sensible to this renormalization scale, is violating this statement 
even for Eiab > 20 MeV. This result implies that some essential physical information for 
such observable is being lost when using /i ~ rup, and as we observe for other observables, 
there is no danger in increasing the value of fi. Another important observation, related to 
a statement of caution made in ref. [7], is that we are not unnecessarily complicating the 
theory, or introducing irrelevant degrees of freedom when considering /i ^ cxd, as it is pos- 
sible to occur in other EFT approaches. Our conclusion for a practical calculation with 
the leading order term is that one should use a momentum renormalization parameter of 
the order of ~ 4 fm~^ (as suggested in [1]) to obtain qualitative results for the low energy 
NN observables. But, to improve the results for some observables, like the triplet mixing 
parameter shown in Fig. 4 and also the ^Si phase shifts, there is no danger in extending 
/i ^ oo, as the previous results are not affected by increasing such value for the scale. 

In conclusion, we have examined the complete renormalization of the leading order term 
of the neutron-proton interaction. We suggest the use of a T-matrix with a subtracted 
propagator, at an energy — /i^, and we have concluded about an overall insensitivity of 
the low-energy singlet and triplet phase-shifts for momentum scale /i larger than 4 fm~^. 
Such result supports the qualitative suggestion presented in Ref. [1]. However, the re- 
sults for the mixing parameter (see Fig. 4) show that this observable is more sensible to 
the renormalization scale, and therefore more appropriate than the diagonal phase-shift 
results to be described by an EFT. Our results for this observable show that fi —>■ oo is 
necessary to have a better agreement with data. And notice that a larger value of fi is 
also consistent with the previous results obtained for the phase-shifts. We can conclude 
that, in our renormalization scheme, fi vanishes as a physical parameter, and the theory 
is renormalization group invariant in the limit /x — * oo [^. 



By excluding the non-diagonal matrix element one could reach a wrong conclusion about ji. 
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It has been achieved with this work that the renormahzed T-matrix for the leading or- 
der interaction has a well defined limit for the physical observables for /z — >■ cx3 with the 
constraint that T{—^'^) being equal to OPEP plus contact interaction, which also is sup- 
ported by the data in the triplet channel. The following step to be investigated is the 
inclusion of the next order term in the EFT expansion of the NN interaction, i.e., the 
Laplacian of the delta potential. In this case, more subtractions are required to construct 
the renormahzed T-matrix, starting from the T-matrix (for the leading order) obtained 
in this work. 

We thank Drs. A. Delfino, M. Malheiro and P. Sauer for useful discussions. This work 
was partially supported by Fundagao de Amparo a Pesquisa do Estado de Sao Paulo 
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Fig. 1. Neutron-Proton phase shifts for the ^Sq state, 6o^s- The sohd hne is the model result, for 
different parameters fi; the dotted line gives the fitting to data from Ref. [15]. 
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Fig. 2. Neutron-Proton phase shifts for the "^Si state, (5o,t- Calculations are presented for ^u 
10 and 30 fm~^; the dotted line gives the fitting to data from Ref. [15]. 
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Fig. 3. Neutron-Proton phase shifts for ^Di state, S2,t- The sohd hne is the present calculation, 
for different parameters fi >4 fm~^; the dotted line gives the fitting to data from Ref. [15]. 
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Fig. 4. Results for the mixing parameter (in degrees) given by eq. (27), for /x 
fm"^. The dotted hne gives the fitting to data from Ref. [15]. 
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